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Abstract 

I present a weak lensing sensitivity estimate for upcoming high redshift (epoch of 
reionization and beyond) 21cm surveys. Instruments such as PAST, LOFAR and 
SKA should be able to measure the weak lensing power spectrum to precisions far 
exceeding conceivable optical surveys. Three types of sources are detectable, which 
include the re-ionization stromgren spheres, large scale structure, and minihalos. 

Tomography allows the measurement of the time evolution of the dark matter 
power spectrum. Raw sensitivities allow measurement of many cosmological param- 
eters, including dark energy, neutrino mass and cosmic equations of state, to percent 
accuracy. It also has the potential for inflationary gravity wave measurements. 

Ultimate limits may be achievable through radio observations of 10 18 minihalos. 
Inflationary Hubble parameters Hi down to 10 _9 Mpi anc k can be detected through 
this effect. Second order effects may also be observable, allowing tests for backreac- 
tion and the quantum mechanical origin of perturbations. 

Key words: Cosmology-theory-simulation-observation: gravitational lensing, dark 
matter, large scale structure 
PACS: 98.80.Es,95.35.+d 



1 Introduction 



Weak gravitational lensing has recently been established as a clean tool to 
probe the distribution of all matter using only its gravitational effects. This 
can also be modelled from first principles, and the matter p ower spectrum has 



already been measured to some accuracy ( Pen et al. . 2003). 



Current measurements of the cosmic microwave background have resulted in 
power spctra that are accurate at the percent level when appropriately binned. 
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This constrains some combinations of cosmological parameters to similar per- 
cent level accuracies. But there are many degeneracies in standard cosmo- 
logical parameters arising from purely CMB data. Using the WMAP data 
alone, even unpo pular regimes are al lowed, such as a closed model with low 
Hubble constant ( Spergel et al. . 2003| ). The authors dismiss such models as 
"unreasonable" , primarily due to the non-standard required value of Ho = 35. 
Nevertheless, one can find references in the rece nt literature which predicted 



such low Hubble constants! Bartlett et al. . 1995| ). 



Some parameters, such as the absolute baryon density VL^h 2 are very well 
determined by the CMB alone, while many other physical parameters require 
external data to constrain accurately, and thus have much larger uncertainties. 
Some such examples include the matter density of the univere, equation of 
state, or neutrino mass. Precision low redshift measurements are also needed 
to break these degeneracies. Galaxy distributions have been used successfully 



as a calibrator ( Spergel et al. . 2003| ). but its precision may have fundamental 



limitations. The external calibrators tend to be limited by systematic errors 
in the models, whose error bars are difficult to quantify. 

Robust low redshift degeneracy breaking measurem ents have already come 
from weak gravitational lensing ( Contaldi et al. . 20031 ). Future lensing surveys 
will signficantly improve on the statistical accuracy, which is currently at the 
10% level. Accuracies greater than 1% are in progress ((Van Waerbeke and Mellier 
20031 ). Much of the most interest ing information, such as the equatio n of state 
of the universe or neutrino mass (|CooravLlT999llFukugita et all [2000). requires 
even more precise measurements. Multiple redshift measurements would be 
invaluable. Large physical scales are important to measure the linear power 
spectrum, which are not affected by baryon feedback. All of these are very chal- 
lenging for the current weak lensing strategy. The surface densities of galaxies 
on the sky is limited, and large sky coverage is expensive. 



The ideal source for weak lensing maps should be at high redshift, and have 
structure on small scales. The CMB satisfies the first criterion, but not the 
second. Faint galaxies partially satisfy both. The epoch of reionization uni- 
verse is a natural candidate to use as a source screen to measure gravita- 
tional lensing. It is far away, emits brightly in the hydrogen hyperfine tran- 
sition line, has structures on many scales ranging from several arcminutes to 
less than a milliarcsecond, a high surface density, and redshift information 
on each source. Several experiments are under construction, in cluding PAST 
(http : / / astrophysics . ph ys . emu . eduT^ jbp /past 6.pdf), LO FAR ( Kassim et al 



2000) and SKA is in the design phase ((van Haarleml . 120 00) . These instruments 
will measure this source population to high accuracy. Several more experi- 
ments, including T-REX (http:/ /orion. physics. utoronto.ca /sasa/Download/poster/ 
casca_poster.pdf) and CATWALK (ftp:/ /ftp. astro. unm. edu/pub/users/john/AONov03. ppt ) 
may result in earlier detections of re-ionization effects, but at lower signal-to- 
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noise. The Canadian Large Adaptive Reflector CLAR ( Carlson et al. . 2000) 



will also measure this redshift range at lower spatial resolution. In this paper 
I present sensitivity estimates for these experiments, and continue to estimate 
the cosmological limits if one exhausted the information. 

Some of these instruments do not have finalized designs and the cosmic reion- 
ization history is not well known. Therefore we only make rough estimates, 
which we expect are probably good to a factor of two. 



2 Pre-reionization Weak Lensing 



We first estimate the strength of gravitational lensing f or sources at r edshift 
z s = 9. We use the WMAP cosmological parameters ( Spergel et al. . 2003|) 



with Q = 0.27, Q\ = 0.73, a$ = 0.84. The angular diameter distance to 
redshift 9 is 2.23c/ Hq = 6.7/i _1 comoving Gpc. One arcsecond is 33/i _1 kpc 
(comoving). 

Figure 1 shows the lensing convergenc e (n) power spectrum com puted by 
integrating the Limber equation over the lPeacock and Doddsj Jl996) nonlinear 



matter power spectrum. The horizontal axis represents the angular scale, with 
9 ~ 2tt/1. The vertical scale is the variance in the k power, which corresponds 
to the magnification factor of point sources in the weak lensing limit. The plot 
also shows the Poisson noise from PAST/LOFAR and the second generation 
PAST+/SKA. 

The most prominent feature is the peak at / ~ 10 4 . This corresponds to an 
angular scale of two arc minutes. Current lensing surveys have signal to noise 
less than unity at the peak of the spectrum, which prevents the construction of 
direct dark matter maps. At these higher redshifts, the lensing power is a few 
parts in a thousand. The rms amplification or shear is several percent, which 
makes its measurement more tractable than the one percent level fluctuations 
in current cosmic shear surveys. The variations are not linear to the statistical 
accuracy, and strong lensing corrections can be conside red. These have been 
presented, both in perturbative and non-linear form, bv iPen ()2000h . 



Before reionization, most of the baryonic matter in the universe is in the form 
of neutral hydrogen. After z < 30, the majority of that gas gravitationally col- 
lapses with the dark matter into virialized 'mini-halos'. The cosmic microwave 
background has a temperature T CMB = 2.73(1 + z). The gas is collisionally 
heated to the virial equilibrium temperature of the minihalos, which is be- 
tween 10 2 - 10 4 K, so typically hotter then the CMB. The average density of a 
minihalo is 200 times the mean density. At these temperatures and densities, 
the spin temperature is tightly coupled to the kinetic temperature, and decou- 
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Fig. 1. The lensing power spectrum for sources at z = 9. The solid line is the 
shear and convergence power spectrum. Images in the source plane are distorted by 
the square root of the power, which is several percent. The left dotted line is the 
noise level expected from PAST/LOFAR. The dotted line on the right is the noise 
expected from the second generation PAST+ and the SKA. 



pled from the CMB (jMariau et all Il997l : iTliev et all 120021 bonah . In the limit 
of a spin temperature much higher than the CMB temperature, the average 
surface brightness at mean densi t y in s pontaneous emission on the 21cm line 



is T' 



21cm 



23mK ([Madau et all |1997). These minihalos provide a fluctuat- 



ing source screen with a characteristic scale that can be used to measure the 
gravitational lensing effect of the intervening dark matter. In the subsequent 
sections we will discuss the observability of this effect. 
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3 Measuring the power spectrum 



3.1 HighS/N 



Two regimes exist for measuring the weak lensing shear at high redshift. With 
high angular resolution and signal-to-noise, one expects to resolve individual 
minihalos. The characteristic mass is 1O 5 M , and one expects ~ 10 18 such 
minihalos on the sky. Each virialized halo is approximately round, which can 
measure the local shear with a signal-to-noise of order unity. The characteris- 
tic size of a halo is lkpc, corresponding to 30 mas. This requires a baseline of 
D = 13000km, which is resolvable with earth sized baselines. The correspond- 
ing spherical harmonic number is / ~ 10 7 . The temperature of such a halo is 
~ 10 3 K, with a velocity width of about 3km/sec. Virialized halos are about 
a factor of 200 overdense in 3-D. Along the line of sight the radial velocity 
dispersion is comparable to the size of the infall region, so one only expects a 
two dimensional density enhancement, which would be ~ 40. The brightness 
temperature on a halo would then be ~ 0.1K. 

An minimal observing strategy would aim to achieve a signal-to-noise of 1 
on the typical minihalo. A lower signal to noise cannot identify individual 
halos. Higher signal-to-noise is expensive to achieve. The actual accuracy on 
the lensing map reconstruction is boosted by the number of halos along the 
same line of sight that one averages over. We use a simple model to deter- 
mine signal-to-noise. A filled aperture telescope of diameter D corresponding 
to the maximal baseline is used as reference. If the elements are distributed 
uniformly within this diameter, we can treat the actual sensitivity as just a 
dilution corresponding to a telescope with a very small aperture efficiency 
t]a = 4A e fr/("7r.D 2 ) where A c g is the collective effective area of all the elements. 
The noise in the map is then T sys / (rjAVtAu), where we take the integration 
time t to be one year, and the bandwidth the thermal width of each minihalo 
of 3km/sec, Az/ ~ 1.4khz. To obtain a signal-to-noise of unity on each halo in 
a year at a system temperature of 200K, an effective aperture of 200,000 km 2 
is required. As an interferometer, 200 elements with 40 km apertures would 
suffice. A half wave dipole has an effective area of A cS = 3A 2 /(87r) ~ 0.5m 2 , 
so 400,000km of wire are needed to achieve this aperture. Copper has a re- 
sistivity of 1.7 x 10~ 8 f2 ■ m, so a diameter of 0.1mm results in a resistance 
much less than the impedance. The wire requires 40 tons of copper. At the 
date of writing, one pound of copper costs US$1.1, which is a neglible cost 
factor. The amplifiers, interconnects and processing would dominate, but this 
cost reduces by Moore's law. 
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3.2 Low S/N 



Initial experiments such as PAST and LOFAR will measure the distribution of 
gas on linear scales, where the only observable are Gaussian fluctuations, and 
reionization Stromgren spheres. The simplest way to measure the fluctuations 
is to construct a three dimensional map, with the redshift as the third dimen- 
sion. On each point of the angular map, we compute the variance of the density 
field along the z dimension. One way of achieving that is to Fourier transform 
the spectrum, and measure the power in each Fourier mode. In practice, some 
modelling would be needed to achieve appropriate noise weights. There are 
different ways of reconstructing a lensing map from a 'fuzzy' noisy image. If 
only one two dimensional image were known, for example in the CMB, this 
can only be done if we have knowledge of the statistics of the source screen. 
The CMB, by being Gaussian, allows such reconstruction using the connected 
four point functi on which is zero fo r a sta tionary Gaussian process, but is in- 
duced by lensing ( Okamoto and Hu . 20031 ). If the unlensed properties are not 
known, it it presumabl y very d i fficult to disentangle lensing from the intrinsic 
statistics, but also see ICooravl (J200J). 



Our reconstruction will use the intrinsic three dimensional information that is 
available from redshifted 21cm structures. This gains us two kinds of observ- 
ables: the intrinsic correlations of matter are presumably statistically isotropic. 
Lensing will change that anisotropically. Sources that are distant along the line 
of sight but close in projection pass through the same gravitational lensing 
screen, and will have induced cross correlations. It is this latter effect that we 
exploit. We observe the structures at a constant angular scale. Weak lensing 
changes the mapping of angular to physical scale slightly Properties of the 
source generically depend on the physical scales at which they are measured, 
and thus their statistics will also experience a slight change. This is true for 
all objects along the same line of sight, even the distant ones. The variance in 
the beam, for example, will generally increase on small scales. This occurs in 
a correlated fashion for objects over a range in source redshift. The procedure 
now is to take a point in the map, smooth the data along the redshift space 
axis, square that, and compute the mean variance. Now we can ask what the 
optimal smoothing window would be. The smoothing procedure is a convolu- 
tion, which is the same as Fourier transforming the line, squaring each Fourier 
mode, and summing them weighted by a window. Each of the Fourier modes 
is uncorrelated because of stationarity (independent of Gaussianity), so we do 
not need to worry about covariances between fourier modes. We presumably 
want the window which maximizes sensitivity to gravitational lensing in the 
presence of noise. If the data is noise dominated, we weight each point by the 
ratio of signal to noise variance. One needs to know the expected variance as 
a function of redshift, which one can measure from the angular maps. 
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Fig. 2. Matter power spectrum measurement. Solid line: the matter power spectrum 
at z = 9. One expects the 21cm emission to trace the total mass. Dot-dashed line: 
power spectrum of patchy re-ionization with a bias 6 = 4. Dashed lines: noise level 
expected for PAST and second generation PAST+. Dotted line: fractional accuracy 
of power spectrum measurement per logarithmic I bin. 



This results in a two dimensional map a 2 (8 x ,9 y ) whose entries are projected 
variances in the density of neutral hydrogen. We denote the three dimensional 
matter power spectrum by A 2 (k) = k 3 P(k)/(Air). It is shown in figure 2. A 
pixel size 9 maps to a spherical harmonic number I ~ 2n/8. One can think of 
the variance as a measure of A 2 (7) shown in figure 2. Patchy re-ionization can 
significantly boost the 21cm power by modulating the neutral fraction. Figure 
2 shows a model of the boosted power, using a b = 4 bias model on large scales 
ISantoset abi liii. Since the neutral fraction cannot be modulated by more 
than unity, we limited the biased power to the greater of unity and the actual 
matter power. 
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Gravitational lensing has two effects on the variance. It changes the area of 
each pixel in the source plane. Each pixel is a fixed solid angle on the sky. 
The variance measures the power at the pixel scale. As the physical size cor- 
responding to the pixel is increased or decreased by gravitational lensing, the 
same pixel probes a larger or smaller angular scale. From Figure 2 we see that 
moving to larger scales (smaller I) decreases the variance in each pixel. 

This variance map is related to the strength of lensing convergence. The frac- 
tional change in the variance by a fractional change in area is half the loga- 
rithmic derivative of the power spectrum: 

<7 2 (x) K(x)dlQgA 2 (fc) 

(a 2 ) 2 dlogk ' 1 ' 

We expect to evaluate the derivative at the physical scale k corresponding 
to angular scale I of the pixel on the sky. From Figure 1 we see that the 
logarithmic derivative is of order unity. 

The other observable from gravitational lensing is the gravitational shear. If 
one takes a collection of pixels on the sky, one can construct a three dimen- 
sional two point correlation function. Projecting the correlation function along 
the radial direction, we can examine the anisotropy of the two point correla- 
tion function. In the absence of lensing, the 2-D two-point correlation function 
is statistically rotation invariant on the sky. Isovariance lines of the correlation 
function are circles. Gravitational shearing will distort the correlation function 
on each patch to appear anisotropic. The ellipticity is a measure of the gravi- 
tational reduced shear. The signal-to-noise of the two methods is comparable. 
The benefit of the shear procedure is that no prior knowledge of the power 
spectrum is required. 



3.3 LOFAR/PAST 



We examine a rough sensitivity estimate for PAST and an idealized LOFAR 
configuration. For this purpose, we measure the gas power spectrum at z ~ 
9. We assume a system temperature T sys = 400K which is dominated by 
the galaxy, and a baseline of 2km as is planned for PAST. The 21cm line 
becomes 2.1m, and the telescope has an effective resolution of 1/1000 radians 
if it were a single dish. We can approximate this to be a spherical harmonic 
number of 3000. The angular scale corresponds to three arcminutes, and 7 
h^ 1 Mpc c qmoy i ng. T he redshift difference for the same radial separation is 
Az = 0.04 ([Penl Il999h . which is a bandwidth of 500 khz. Initial PAST will 



have an effective aperture of 40,000 m 2 (140Mhz/i/) 2 , which corresponds to 
an aperture dilution tja of about 1/100. This results in a noise level after a 
year's operation of T sys / (r/AVtAv) ~ lOmK per synthesized beam of three arc 
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minutes. Current models place the size of ionization spheres at similar scales 
as this resolution, so one expect of order unity signal to noise. Radially, the 
resolution is much better, so one could realistically expect 100 independent 
radial modes in each synthesized beam across the epoch of reionization. If 
each radial mode represents a signal to noise of order unity, our sensitivity 
is about Ak ~ 0.1. The square is the lensing noise shown in Figure 1 at the 
corresponding angular scale / ~ 3000, and scales as white noise. 

For gravitational lensing, one optimizes throughput by using the scale where 
the signal-to-noise is of order unity. Our procedure is to square the density 
map at that scale, and use the changes in this variance (square) to estimate the 
lensing strength k. The power spectrum of K is now the two point correlation of 
our density variance, which is related to a four point function of the underlying 
density field. Because of our large radial average of squares, we assume that 
the contribution from the intrinsic spatial connected four point function is 
neglibible. With this procedure our lensing noise power scales oc I 2 . 



The actual sensitivity of the telescope to neutral hydrogen density structures 
as a function of scale in the source plane is a different calculation. It depends 
on the actual layout of the elements, and the number of elements at a given 
separation. If we keep the collecting area, and just shrink the baselines, the 
aperture dilution rj^ decreases as the square of the separation, while the radial 
scale is proportionate to the angular scale. The noise scales as # _2 ' 5 . At dou- 
ble the angular scale, corresponding to six arc minutes, this corresponds to a 



noise level of 13 mK J week/t, so PAST can map out larger Stromgren spheres 
directly. It has been proposed that the ionization s pheres up to 50 Mpc hav e 
been observed in existing QSO absorption line data (jWvithe and Loe b. 2004). 
These would be observable at high signal to noise with PAST. LOFAR has a 
similar effective aperture, and we would expect a similar sensitivity. The sec- 
ond generation PAST (which we call PAST+ in this paper) is expected to have 
ten times the resolution and one hundred times the effective aperture. This is 
four times larger than the SKA, so for estimation purpose we clump these two 
instruments into the same class, as shown by the more sensitive noise lines in 
Figure 1. To estimate the total sensitivity on the full field power spectrum, we 
assumed that the PAST array instantaneously measures 100 square degrees of 
the sky, and we measure a 30% redshift depth. For PAST+/SKA we used 10 
square degrees. 



The problem can be broken into two parts: one must measure the power spec- 
trum of the 21cm emission, and then look for spatial variations in the power. 
Figure 2 shows the expected matter power spectrum. We show the noise of 
PAST/LOFAR as the dashed line. We also plotted the statistical accuracy of 
the power spectrum measurement as dotted lines. The matter power spectrum 
itself can of course also provide cosmological information. Unfortunately, this 
will again be limited by the complication of gastrophysical processes such as 
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radiative transfer and star formation, which will modulate the spatial distri- 
bution of neutral hydrogen. 



4 Processing Challenges 



Several problems a rise in the endeavor to map gravitational lensing to such 
an ambitious scale. Oh and Mack ( 2003h have shown that fluctuations in syn- 
chrotron emission from ionized gases can outshine the fl uctuations from the 21 
cm emission. As pointed out bv iFurlanetto et al.l (J200J), this can probably be 
addressed by removing all power-law spectrum spatial fluctuations. A similar 
requirement holds for the modelling of radio point sources. These sources are 
valuable to calibrate the point-spread function. A promising strategy is to use 
the brightest point sources in the map to calibrate the spectrum. The radio 
sources are emit synchrotron radiation from relativistic electrons, which has 
no spectral structure. The task of looking for spectral features now consists 
of comparing the spectrum of points on the map with that of nearby bright 
sources. This can be done to high accuracy. 



Computationally, a brute force correlation of a very large number of dipoles 
can be expensive. A direct correlation for a full sky synthesis with N dipoles 
with involves 0(N 2 ) correlations. Some designs allow beam forming within 
each synthesis node. For nodes with n dipoles placed on regularly spaced pat- 
terns, the cost to form the n independent beams can be reduced to O(nlogn), 
so the total correlation cost is 0((N 2 /n) logn). Measuring the all-sky actual 
power spectra can also be computationally challen ging, especi ally if correc- 
tions to the weak lensing approximation is desired (IPenl. l2000h. The general 



theoretical framework is well understood ( Okamoto and Hul . 120031 ) . but a de- 
faced computational pipeline will require significant software development 
f|Pe^l . |200? 



Since the array is sparse, not all projected baselines are measured simultane- 
ously. Earth rotation fills in the map. To perform statistics on 10 18 sources of 
course will require a correspondingly large catalog. The information is mea- 
sured in the exabytes. With today's hard drives, this would cost billions of 
dollars. But storage and computation both follow Moore's law, and will be- 
come affordable in the forseeable future. Certainly many challenges must be 
overcome to implement the ultimate array. 
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5 Cosmological Constraints 



We discuss some of the cosmological constraints achievable with pre-reionization 
21cm emission. The cosmic shear can be measured to unity signal-to-noise for 
I < 500 for PAST/LOFAR, I < 10000 using PAST+/SKA and I < 10 7 using 
the ultimate detector. 

Nominally, this results in an accuracy of projected dark matter power spec- 
trum determinations to better than 10~ 4 for PAST+/SKA, and better than 
10~ 7 in principle. If the only unc e rtaint y came from the neutrino mass, we can 
apply equation (6) from Coorav ( 1999j ) and derive a nominal accuracy on the 



neutrino mass of O.lmeV for PAST+/SKA and \i eV as the ultimate limit. 
This is already significantly more accurate than the mass differences between 
neutrinos, and it should be possible to measure the individual masses of each 
of the neutrino generations separately. These masses can then be compared to 
experiment. 

Over a fiducial range of source redshifts 8 < z < 10, the angular diameter 
distance varies by 5%. For a lens half way between us and the source z ~ 1.6, 
the lensing strength changes by 3%, which is easily measurable. Similarly, to- 
mography allows the reconstruction of the power spectrum at two different 
redshifts. For a lens screen 1/3 and 2/3 to the source screen (redshifts 0.9 and 
2.7 respectively), the relative amplitudes are only 90% degenerate, loosing 
only a factor of ten in accuracy for determining the difference in amplitude. 
The change in power spectrum can be measured to an accuracy of 10~ 3 , which 
allows very accu rate dynamical constraints on the equation of state of the uni- 
verse f[Hull2002h . Most of the proposed cosmological constraints by precision 



measurement of the dark matter power spectrum can be implemented in this 
precision tomographic surve y. Baryon oscillations have be en suggested as an 
example of a standard ruler ( Blake and Glazebrook , 2003| ). 



Inflation generically predicts the existence of gravity waves produced by Hawk- 
ing radiation from the de Sitter horizon at the time of inflation. The ampli- 
tude at the quadrupole is C\ ~ (Hj/M p \ anck ) 2 , where Hj is the Hubble pa- 
rameter during inflation. Using I 3 ~ 10 12 effectively independent sources i n 
PAST+/SKA, we can constrain Hi < 10" 6 M p i anck (|Dodelson et all l2003h . 



The 21cm fluctuations should be at least 2 orders of magn itude more sensi- 
tive than the galaxy statistics used in lDodelson et al. (|2003|) . For the ultimate 



survey, we expect to see 10 18 minihalos on the sky. If one resolves each, and 
measures weak lensing, we can measure inflationary expansions all the way to 
a Hubble parameter at inflation of Hi ~ 10~ 9 M, 



planck • 



Such a precise measurement of the power spectrum would also unveal second 
order effects. Perturbations at inflation have amplitudes of order O(10~ 5 ). Var- 
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ious second order effects enter at an amplitude of 1CT 10 , and any measurement 
of the power spectrum more accurate than one part in 10 5 will depend on the 
second order initial conditions. 

These estimates are all idealized, and signficant work is needed both in the 
instrument design and the processing stages to achieve these sensitivities. 



6 Conclusions 

I have presented sensitivity estimates for epoch of reionization 21cm weak 
lensing surveys. We find that a precise mapping of the projected dark matter 
structure on the sky can be achieved with several projects under construction 
and planning, including PAST, its second generation PAST+, LOFAR, and 
SKA. In principle, it is also possible to image up to 10 18 minihalos with an 
extremely large dipole array. Such a future step would require an effective 
aperture of > 10 5 km 2 . The overwhelming statistical accuracies would allow 
accurate measurement or stringent upper bounds on primordial gravity waves, 
and second order effects from the epoch of inflation. 

Even with the upcoming experiments, exquisite maps of dark matter and 
its evolution is possible, at much higher precision than any other planned 
or proposed procedure. The two dimensional dark matter distribution can be 
measured to a signal-to-noise of better than unity at scales down to I < 10000. 
With the leverage of a redshift range 8 < z < 10, the evolution of the power 
spectrum, dark energy and neutrino mass can be measured to better than 
percent accuracies. 

I thank Martin White for helpful discussions and Pengjie Zhang for the lensing 
code. 
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